Generation of Large Number-Path Entanglement Using Linear Optics and 
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We show how an idealised measurement procedure can condense photons from two modes into 
one, and how, by feeding forward the results of the measurement, it is possible to generate efficiently 
superpositions of components for which only one mode is populated, commonly called "NOON states" . 
For the basic procedure, sources of number states leak onto a beam splitter, and the output ports 
are monitored by photodetectors. We find that detecting a fixed fraction of the input at one output 
port suffices to direct the remainder to the same port with high probability, however large the initial 
state. When instead photons are detected at both ports, Schrodinger cat states are produced. We 
describe a circuit for making the components of such a state orthogonal, and another for subsequent 
conversion to a NOON state. Our approach scales exponentially better than existing proposals. 
Important applications include quantum imaging and metrology. 



The fundamental limits to optical detection for metrol- 
ogy and imaging are quantum mechanical [l[ • Of partic- 
ular interest for reaching such quantum limits are path- 
entangled states of photons of the form \N0) + e^\QN), 
in a basis of photon-number states, commonly referred to 
as "7V00-/V" states. A variety of applications have been 
suggested @. For lithography [l| and microscopy 
NOON state light would be used together with multi- 
photon absorbers to achieve enhanced resolution. This is 
because the de Broglie wavelength for an AT-photon state 
is a factor 1/N smaller than the wavelength associated 
with the single photon, and the absorption rate scales 
linearly with the incident intensity, rather than as the 
N th power. Regarding applications to precision metrol- 
ogy, whereby an interferometric setup is used to mea- 
sure small phase shifts, NOON states achieve the Heisen- 
berg limit, for which the phase uncertainty scales as 1/N 
[H 13 1 and entanglement is a fundamental requirement 
for achieving this limit. It has been rigorously demon- 
strated that the cost of improving sensitivity (without 
using entanglement) is higher intensities or longer coher- 
ence times [8]. Classically the shot-noise limit applies, 
attained for example by laser light, for which the uncer- 
tainty scales as 1/yN, already a restraint in applications 
such as magnetometry [9( and gyroscopy [10| . 

However, building a source of NOON states beyond two 
photons is challenging. Three, four and six photon ex- 
periments have been reported [ll|, E 03 , but only in the 
first two references were NOON states generated. In the- 
ory a source could be made using a nonlinear crystal [l4| . 
However, the required optical nonlinearity is not readily 
available. An alternative is a non-deterministic approach 
using linear optics, wherein the desired state is generated 
on condition of a specific outcome at photodetectors. A 
variety of schemes have been suggested which typicalfy 
rely on conditional destructive interference (l5l . y^l Qj 
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However, so far none of these scales efficiently, that is 
they all share the feature that exponentially decreasing 
success probabilities outweigh the possible gains. Noting 
that quantum algorithms, exhibiting polynomial and ex- 
ponential speedups over their classical counterparts, may 
be implemented scalably in a linear-optics approach [18| , 
we expect that it should be possible to do better. In this 
Article we address this challenge by adapting a concep- 
tually simple measurement procedure. Our method is as 
follows. First, we aim to minimise the negative effects of 
back-action in a sequence of detections, whereby earlier 
measurements affect the outcomes of later ones. Next, 
we engineer output states that closely approximate the 
ideal case. Finally, we exploit feed- for ward, for which cir- 
cuits are actively switched in response to previous pho- 
todetections. Feed-forward is an essential ingredient of 
linear-optics based quantum computing, but is not used 
in previous proposals for engineering NOON states. 

We begin by considering the thought experiment de- 
picted in Fig. 1(a). Here, cavity modes labeled A and 
B are assumed to start with a well-defined photon num- 
ber N. They are each coupled to an external mode by 
a weakly transmissive mirror, and these modes are com- 
bined at a 50:50 beam splitter, and then subject to par- 
tial photodetection. The beam splitter acts to make the 
origin of the photons indistinguishable. When a photon 
is registered at the left or right photodetector labeled 
Dl or Dr, the transformation is given by the Kraus 
operators L = (a — 6)/v2 or R = (a + b)/y/2 respec- 
tively (where a and b are the annihilation operators for 
modes A and B). To obtain the corresponding prob- 
abilities it is necessary to normalise by the total pho- 
ton number prior to detection. We suppose now that a 
string of detections occur only at D r , say (by adjusting 
the path length difference of the cavities between detec- 
tions, with a phase shifter, the same state for the cav- 
ity modes can be obtained in every case). After this, 
the detectors are removed and the system evolves to a 
final state with all the remaining population at the out- 
put ports. Denoting by \iPab) the state of modes A 
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FIG. 1: Two measurement-based procedures, inducing rela- 
tive phase correlations between the principal modes, each one 
with N photons at the start. In (a), population leaks from 
cavity modes A and B into external modes, which are com- 
bined at a beam splitter. Photons are detected one at a time 
by photodetectors Di and D r . In (b), all modes are propagat- 
ing. Beam splitters couple some fraction / from the principal 
modes one and two into ancillae modes three and four, which 
are initially the vacuum. The ancillae are subsequently com- 
bined at a beam splitter, and subjected to number-resolving 
photodetection at D; and D r . 



and B after r initial detections, we find that \iPab) = 
R r \N)\N) /J(N\ (N\ (J#Y (ttj r \N) \N), normalising 

to unity. The probability P C ond of finding all the remain- 
ing photons "condensed" at the right output port (and 
none at the left output port) is as follows, 
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where S = 2N — r denotes the total remaining photon 
number, C denotes a binomial coefficient, and we assume 
that r < N . Evaluating the value of P con d numerically 
for initial states of increasing size, we find that its value 
is determined asymptotically by the proportion of the 
input that is measured. For example, setting r either as 
one quarter or one third of 2N suffices for P con d > 0.6 or 
-Pcond > 0.7, respectively. 

We have found for our thought experiment, that later 
detections tend to strongly reinforce earlier ones. Hence, 
the effect of measurement back-action here is useful for 
state engineering, and in what follows we adapt the mea- 
surement process for NOON state generation. We di- 
vide our analysis into three stages. First, we translate 
our thought experiment into a mathematically equivalent 
procedure based purely on linear optics, and consider the 
general case for which photons are detected at both pho- 
todetectors. The localisation phenomena resulting from 
this measurement process, have been studied extensively 
in the context of the debate over the existence of abso- 
lute optical coherence in common quantum-optical exper- 
iments 19]. It has been demonstrated that well-defined 
correlations in the relative optical phase evolve, for the 
remaining population, and play a central role in the ongo- 



ing dynamics. Hence, in the second stage of our analysis, 
we investigate simple procedures for manipulating phase 
correlations, and relate states with well-defined correla- 
tions in the relative phase and NOON states. Finally, we 
identify a method based on feed-forward to enable NOON 
states to be generated efficiently. 

First, we translate our thought experiment into a 
mathematically-equivalent procedure, based on linear op- 
tics, as depicted in Fig. 1(b). We label this optical cir- 
cuit, Circuit I. Here all modes are propagating, and a 
source is assumed to supply dual Fock states |iV) |iV) to 
the principal modes one and two. Beam splitters of re- 
flectance / couple modes one and two to ancillae modes 
three and four, which are combined at a 50:50 beam 
splitter. They are then measured by number-resolving 
photodetectors labeled Di and D r , where on average a 
fraction / of the input photons are registered. We now 
consider the state \ipi,r) generated in modes one and two 
after I photons are registered at Di and r at D r . Fol- 
lowing Ref. [2(|, it is convenient to adopt a representa- 
tion in terms of coherent states, which are of the form, 
| a) = ||a| exp(i#)) oc J2kLo \/\<x\ 2k /k\ exp (ikO) \k) in a 
basis of Fock states. It has been shown that, 

/•27T P 7T 
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{6 1 + 82) /2 and A = 62 — 8±. The superposition phase 
a takes the value lir, and hence the measurement record 
must be known exactly. The scalar function G(X) is 
given to good approximation by the Gaussian expression 
exp [—{I + r)X 2 /4] . The total photon number is equal 
to S . There are well-defined correlations in the relative- 
optical phase parameter A at values plus and minus Ao, 
determined only by the ratio of / to r. These correla- 
tions are multi-valued whenever photons are registered 
at both photodetectors, and a Schrodinger cat state is 
generated. We can see that cats are generated as a re- 
sult of the symmetry of the setup. Specifically, L and 
R are invariant under an exchange of the labeling of the 
modes, a transformation which reverses the sign of the 
relative phase. The generation of cat states therefore 
also requires precise phase stability between the modes. 
Turning to the source, we see that the state of the input 
can be a mixture of the form ^2 N P^\N)\N) (N\ (N\, since 
Ao is independent of N, and standard linear optical ele- 
ments obey a superselection rule for the photon number. 
Several two-mode squeezing processes strongly suppress 
relative number fluctuations, and hence might serve as 
practical sources of light described by these mixed states. 

For the second stage in our analysis, we identify the 
outputs of Circuit I as examples of quantum reference 
frames — reference frames for a classically defined pa- 
rameter composed of finite quantum resources. Quantum 
reference frames are subject to depletion and degradation 
as they are used, and are currently of interest for proto- 
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cols in the field of quantum information, in which they 
are regarded as a resource [21] . By making an analogy 
to classical phase references we can now identify simple 
ways in which states of the form Eq. |T]) can be manip- 
ulated. For the current purposes we can assume that 
a large number of detections have been performed and 
define, 

IV-oo (A )) oc / <Wexp(-iS6)\a)\aexp(iA )), (2) 
Jo 

where a — \a\e l6 , for a state with a total photon num- 
ber S and a relative phase of A (assumed to be nor- 
malized) . Relative phase correlations between more than 
two modes are transitive and are transformed additively 
by phase shifters. A phase reference can be extended to 
additional modes by combining it with the vacuum at 
a beam splitter. As an example, a 50:50 beam split- 
ter, which we denote here by Uh s , beating light in a 
Fock state with S photons against the vacuum yields, 
U hs \S)\0) oc |Voo(0)), and U hs \0)\S) cx There- 
fore, we see that a simple circuit, consisting only of a 
beam splitter and a phase shifter, can convert a cat state 
generated by Circuit I to a NOON state, whenever the rel- 
ative phase correlations differ by tt. This happens when 
I = r and Ao = tt/2. We label this circuit, Circuit III 
(anticipating an intermediate process modifying the cat 
states for the general case). 

Before proceeding to the final stage of our analy- 
sis, we consider a simple AOOiV-state generator, that 
attempts to convert every cat state generated by Cir- 
cuit I using Circuit III. This method might be expected 
to yield close approximations to NOON states, when- 
ever the the relative phases of the cat state are close 
to plus and minus 7r/2. The situation is summarized 
in Fig. 2(a). To measure the quality of the output 
state we adopt the fidelity, denoting it by F. For 
the measurement-induced condensation, considered at 
the start, F takes the same same value as P C ond- For 
schemes generating NOON states, it is necessary to ac- 
count for the phase of the superposition, and we define 
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F 



where 



output/ 



max (<0S| + exp (-10) <S0|) ^output) A 

S is the total photon number of the state \ip, 
Evaluating F for our TVOOTV-state generator, when Cir- 
cuit I generates a cat state with relative phase compo- 
nents at ±Ao and total photon number S, we find to 
first approximation that F ~ cos 2S [(Ao — 7r/2) /2]. As 
with other proposals, this scheme in fact scales exponen- 
tially poorly whenever the relative phase correlations are 
less than tt apart, as is typically the case. Inspecting 
the overlap for different relative phase components, as 

in Eq. Q with total photon number S 1 , we find that 

S 



<lMAiWoo(A 2 )) 



The poor 



cos [(A 2 - Ai) /2] 
scaling can be attributed to the non-orthogonality of the 
cat state components. 

We now proceed to the final stage of our analysis. Our 
previous iVOOiV-state generator is effective when Circuit I 




FIG. 2: (a) and (b) illustrate complete NOON state generators 
in outline. Circuit I produces Schrodinger cat states with two 
relative phase components non-deterministically (green). The 
generators are terminated by a fixed unitary (yellow) circuit 
consisting of beam splitters and phase shifters. In (b) the cat 
states are corrected, using an additional measurement process 
conditioned on the previous detection outcomes. 



generates cat states with components which are orthog- 
onal. However, this occurs with low probability. It is 
not possible to improve the situation with any combi- 
nation of (idealised) beam splitters and phase shifters, 
since these implement unitary transformations. Hence, 
we now devise a circuit, labeled Circuit II, to make input 
cat components orthogonal, using additional processes of 
measurement and feed-forward. This more sophisticated 
scenario is depicted in Fig. 2(b). To identify a suitable 
circuit, we investigate how a beam splitter transforms 
phase references, starting with two classical fields. Here 
the field in each mode is represented by a complex num- 
ber, with the square amplitude corresponding to the in- 
tensity, and the phase to the optical one. A 50:50 beam 
splitter, configured so as not to impart additional phase 
shifts to the modes, outputs two classical fields described 
by the sum and difference of the values for the inputs (al- 
tering both the square amplitudes and the phases). If the 
input has a relative phase of or tt, and equal intensities 
for each mode, the population is transferred entirely into 
one mode. On the other hand, if the input has a relative 
phase of plus or minus 7r/2, and equal intensities in each 
mode, the relative phase and intensities are preserved. 
Moving to the quantum case, we consider the action of 
the beam splitter for a state, defined as in Eq. with 
a relative phase of Ao, a total photon number S, and an 
intensity S/2 in each mode. Computing the final state 
explicitly, we find a scenario similar to the classical case, 

i-2-k 

Uhs\i>oo) oc / d6 exp(—iS9) x 
Jo 



h exp (i0)) | y/I 2 exp [i {9 ± tt/2)]) . (3) 

This has an intensity SI\/ (I x + J 2 ) = [1 - cos(A )] /2 
in mode one and SI2/ (h + I2) = S [1 + cos(Ao)] /2 in 
mode two, and a relative phase of plus 7r/2 when < 
Ao < 7r/2 and of minus tt/2 when — tt/2 < Ao < 
(we consider cases for which the intensity is increased in 
favour of mode two) . 

The symmetry of the beam splitter transformation 
makes it useful for altering the cats generated by Circuit 
I, so that the relative phases are different by tt. How- 
ever, it creates a difference in the intensities between the 
modes. To correct this, we propose beating mode two 



4 



against the vacuum, so as to to move the difference of 
the intensities to an ancillary mode, which can be re- 
moved by a photodetection. This method depends criti- 
cally on feeding forward the result of the detections per- 
formed by Circuit I, so that a variable beam splitter can 
be set according to the value of Ao. A variable beam 
splitter can be implemented with 50:50 beam splitters 
and variable phase shifters. The cost of correction is a 
decrease in the total photon number, which varies non- 
deterministically. As can be seen from Eq. , a fraction 
of cos (Ao) of the photons are lost on average. Overall, 
Circuits I through III constitute a complete AOOA-state 
generator. Additional mathematical analysis is given in 
the supplementary online text. Runs for which Circuit 
I fails to generate a Schrodinger cat state, or too many 
photons are lost in the detection process are discarded. 
The fidelities at the output are, on average, 0.87, 0.94 or 

0. 98, when a fraction of one third, one half, or two third 
respectively of the input photons are detected by Circuit 

1. Higher fidelities are possible when the photon number 
at the input is small. If allowance is made for sufficient 
input photons to be detected by Circuit I, and a further 
half to be detected in Circuit II, the probability of failure 
is not too large. 

Finally, we suggest some possibilities for experimen- 
tal implementation. For the source, we propose an op- 
tical parametric oscillator setup for which the two mode 
squeezed output of an optical parametric amplifier is en- 
hanced by a cavity [22| . Note, however, that the cur- 
rent purposes require twin beams of a much lower inten- 
sity than is typical in many experiments, and that the 
beams must be rendered frequency degenerate. Tech- 
niques of feed-forward and photodetection are being de- 
veloped with a view to quantum information technolo- 
gies [23l . |24| . For the source, an important problem is 
imperfect correlation between the modes. If, for exam- 
ple, two independent lasers of equal intensity provide the 
input, the scheme generates the intended relative phase 
correlations, but no entanglement (25[. Photodetectors 
are subject to loss and dark counts. Losses will act to 
degrade the source, reducing the relative number corre- 
lation and increasing the uncertainty in the total photon 
number. Dark counts are more problematic, mixing over 
the phase for the superposition in Eq. Jl}. An alter- 
native suggestion is using trapped bosonic atoms. One 
possibility might be to work in a regime for which the 
atomic wave-packets are much longer than the typical 
scattering length, as proposed in Ref. [26j]. Another is 
to use Bose-Einstein condensates, for which a variety of 
coherent operations have been demonstrated. Number- 
resolved condensates might be obtained from the Mott 
Insulator phase, while relative-number squeezing can be 
achieved by different techniques. 

In conclusion, we have proposed for the first time a 
linear-optics based scheme that generates large A00A 
states efficiently, the photon number at the output scal- 
ing with that of the source — all the while maintaining 
high fidelities, high success probabilities and a fixed num- 



ber of circuit components. As well as being of immediate 
interest for a range of applications, our results have con- 
nections with other topics. For example, the scaling we 
derive for our measurement-induced condensation pro- 
cedure is of relevance to the study of the interference of 
light from independent sources and localizing relative op- 
tical phase, phenomena with analogs in different physical 
systems [27j j . We have left as an open question the extent 
to which the scaling can be attributed to Bose statistics. 
Regarding our AOOA-state generators, the creation of 
macroscopic entangled states is of interest for exploring 
the quantum-classical transition. Finally, our study of 
Schrodinger cat states may have application to quantum 
computing, where Schrodinger cat states, defined for one 
mode only, have been proposed to encode qubits, which 
may be manipulated using standard experimental tech- 
niques [2^ . 
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Supplementary Material: Methods 



evolves to a cat state according to, 



In these supplementary notes, we provide further anal- 
ysis of our AOOA-state generator, consisting of Circuits 
I, II and III, as depicted in outline in Fig. 2(b). First, 
we specify notation for beam splitters, phase shifters and 
states with well-defined relative phase correlations. For 
the lossless beam splitter, we choose a notation which 
makes explicit the "rotation" performed by such a device. 
A beam splitter with transmittance r and reflectance 
(1 — r) acts to transform the annihilation operators dj 
for modes labeled j, according to the relations, 



cos(7) - sin(7) \ / 01 
sin(7) cos(7) J I 6 2 



with angular parameter 7, where t = cos 2 (7) and < 
7 < tt/2. We denote this transformation (7), and we 
include, where necessary, phase shifts of x at the input 
port and — x at the output port of the first mode, so that 

Uhs(j,X) = exp(jexp(ix)didl-'rexp(-ix)d{6 2 y For 

example, Uhs (arccos [\/t), tt/2) corresponds to a sym- 
metric beam splitter. We denote a phase shift trans- 
formation on mode j, exp (idj ' djx) , by U ps (x)- For a 
state defined, as in Eq. (2), with a total photon number 5* 
and relative phase Ao, it is helpful to incorporate a phase 
factor exp (— iSAo/2) into the normalisation (making the 
definition symmetric between the modes) . We then adopt 
the following notation for a normalised Schrodinger cat 
state, 

l^cat (A , A)) cx l^oo (Ao)) + exp (iA) |Voo (-A )) , 

having components with relative phases plus and minus 
A , a phase for the superposition A (with the overall 
normalisation constant assumed positive). 

Next, we elaborate on the sequence of operations per- 
formed by our AOOA-state generator. We assume the 
final state should have at least P photons, and that the 
correlations in the relative phase are ideal. For the first 
step, Circuit I, depicted in Fig. 1(b), implements the 
transformation, 



\l,r)(l,; 



asm 



asm 



Vf' 



conditioned on the detection of I photons in mode 3 and 
r photons in mode 4. A dual Fock state from the source 



Circuit I :|A) 1 |A) 2 |0) 3 |0) 4 — > IV'cat (A A)) lj2 . 

This cat state has relative phase components with values 
plus and minus Ao = 2 arccos \Jr/(l + r) , and total 
photon number 2N — I — r. Runs for which I or r are 
zero must be discarded. It has been shown that values 
for the relative phase are generated with approximately 
equal frequency across the range [25[, and hence these 
failure events do not affect the scaling of the generator. 

Next, Circuit II acts to transform the relative phase 
correlations, to plus and minus tt/2, in every case. When 
r > I, the relative phase correlations lie in the range 
[—tt/2, tt/2], and a 50:50 beam splitter acting on the prin- 
cipal modes corrects the relative phase correlations, while 
increasing the intensity in mode two (and decreasing it in 
mode one). To achieve the same outcome when I < r, we 
suppose that a phase shift of tt is applied in advance (on 
either mode) . This transforms the cat state generated by 
Circuit I as, 

U ps (tt) |^c a t (A (Z,r),Z7r)) cx \ip cat (A (r,l),nr)) . 

Next, a beam splitter, with transmittance 
[1 — cos (A )] / [1 + cos (Ao)], transfers the differ- 
ence of the intensities to the ancillary mode five. A 
circuit for implementing the variable beam splitter is 
given by the relation, 

U bs (7)2,5 = 

U hs (tt/4, tt/2) 2i5 U ps ( 7 ) 5 U ps (-7)2 ^bs (tt/4, -tt/2) 2 5 . 

A photodetector measures Q photons in mode 5. Overall, 
Circuit II implements the transformation, 

\Q) (Ql 5 Uhs {arccos [tan (A /2)]} 2 5 U hs (tt/A) 12 . 

The cat state evolves as, 

Circuit II : 

|V>cat (Ao, Itt)) — > l^cat [tt/2, ln+(2N-l-r- Q)tt/2\) 

We derived a full probability distribution for the out- 
comes of the photodetection performed by Circuit II, 



Prob(Q = 0, • • • , S - 1) = 
Prob(Q = S) = 



; - J ~s~7~a I S °Q [l-cos(A )] S - g cos<2 (A ) 
1 + (-1) cos 6 (A ) 



1 + (-!)' 



l + (-l)W (Ao) 



■cos s (A ) , 
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where S — 2N — I — r is the total photon number prior 
to detection, and C denotes a binomial coefficient. This 
probability distribution is approximately binomial, and 
the expected number of detections is S , cos(A ). If too 
many photons are lost in Circuits I and II the run must 
be aborted. The probability of this can be made small 
by taking P ~ N(l — /). In principle, excess photons can 
be removed by an additional process, similar to Circuit 
I. 

Finally, Circuit III implements the unitary transfor- 
mation, 

U hs (7r/4,7r) 12 f/ ps (tt/2) 2 . 
The corrected cat state evolves as, 



Circuit III : 

IVcat (tt/2, l7r+(2N-l-r-Q)w/2)) - 



|P,0) + (-l)'|0,P), 



yielding the desired NOON state, with P = 2N-l-r-Q. 
It may be noted that the superposition phase for the 
iVOOiV-state at the output depends on the measurement 
record at the photodetectors. When I < r, the additional 
phase shift in Circuit II causes this phase to be rir rather 
than lit. 

Next, we estimate the fidelities of the states produced 
by our NOON state generator, and clarify its behaviour 
for large photon number. To do this, we first compute 
the fidelity for one component of a cat state generated 
by Circuit I, which we denote by |^g(Ao)). We assume, 
as in Eq. (1), that the function G(X), describing the 
localisation of the relative phase, assumes its Gaussian 
asymptotic form. Note that the rate of localisation is 
faster when phase correlations evolve at more than one 
value. We assume that the state at the input is the dual 
Fock state |AT) 1 |AT) 2 , and that a total of D = I + r de- 
tections have occurred. Then, 



I + r = 5, and range from 0.96 to 1 when I + r = 6, 
while for input state 1 15) 1 1 15)2 the values range from 
0.92 to 0.96 when I + r = 15, and from 0.96 to 0.99 when 
I + r = 20. 

Finally, we performed a complete numerical simulation 
of the iV00-/V-state generator, to verify that Circuits I, II 
and III work together as predicted. In particular, it was 
necessary to check that Circuits II and Circuit III func- 
tion as expected when the relative phase correlations for 
the cat states are not perfectly well-defined. The results 
are shown in Fig. 3. Each point in the plots corresponds 
to a particular choice of input state, and measurement by 
Circuit I. The height corresponds to the expected pho- 
ton number for the output iVOOiV state, and the color 
to its fidelity. Averages are taken over all possible out- 
comes to the third photodetection performed by Circuit 
II. For comparison, the mesh shows the predictions of 
the preceding analysis. Good agreement is seen between 
these analytical predictions and the numerical results. 
However, inspection of individual outcomes in Circuit II 
reveals that the high fidelities are not maintained in ev- 
ery case. Roughly speaking, improbable outcomes were 
often found to have low fidelity. 



(Vg(Ao)IVoo (A )> 



/:/ 2 /2 dAcos*(#)exp(-^ 



f-i% 1% dAdA' coss ( A-Ai) cxp [_d (A ,2 + A2)] 



v [ 2 m-i] 

where S — 2N — D is the total photon number. This 
result was derived assuming that S and D are not small. 
The value for the fidelity depends only on the ratio of 
detections to input photons. For example, when D/2N 
is one half, F — 0.94, and when D/2N is two thirds, 
F = 0.98. To verify this result, we computed numerically 

exact values for the fidelity, (Vvl^cat (Ao(i, r), In)) 

for a range of states \ipi, r ) generated by Circuit I. For 
input state |3)i|3) 2 , the fidelity is 0.94 for (l,r) = (1,2) 
and (2, 1), and anomalously it is 1 for (l,r) = (1, 1). For 
input state 1 5) x 1 5) 2 the values are 0.94 and 0.96 when 



0.88 0.9 0.92 0.94 0.96 0.98 1 

Fidelity 




'0 



FIG. 3: Fidelities (color) and photon number (vertical axis) 
are displayed for outputs of our NOON states generator. Each 
point corresponds to a possible outcome to Circuit I, for which 
D photons are detected. Input states \N)\N)2 are considered 
for N up to 15. Going from left to right, D/2N is one third, 
one half and two thirds. 



